Effects of domain walls on hole motion in the two-dimensional anisotropic t-J model 

at finite temperature 



cn 

O 

o 

(N 



X3 
O 

(N 
> 

O 

O 
cn 
o 



o 



X 



Jose A. Riera 

Instituto de Fisica Rosario, Consejo Nacional de Investigaciones Cienttficas y Tecnicas, y Departamento de Fisica, 
Universidad Nacional de Rosario, Avenida Pellegrini 250, 2000-Rosario, Argentina 
Institute for Materials Research, Tohoku University, Sendai 980-8577, Japan 
(Dated: February 6, 2008) 

The t-J model on the square lattice, close to the t-J^ limit, is studied by quantum Monte Carlo 
techniques at finite temperature and in the underdoped regime. A variant of the Hoshen-Koppelman 
algorithm was implemented to identify the antiferromagnetic domains on each Trotter slice. The 
results show that the model presents at high enough temperature finite antiferromagnetic (AF) 
domains which collapse at lower temperatures into a single ordered AF state. While there are 
domains, holes would tend to preferentially move along the domain walls. In this case, there are 
indications of hole pairing starting at a relatively high temperature. At lower temperatures, when 
the whole system becomes essentially fully AF ordered, at least in finite clusters, holes would likely 
tend to move within phase separated regions. The crossover between both states moves down in 
temperature as doping increases and/or as the off-diagonal exchange increases. The possibility of 
hole motion along AF domain walls at zero temperature in the fully isotropic t-J is discussed. 

PACS numbers: 74.20.-z, 74.25. Dw, 74.25. Ha, 02.70.Uu 



The interplay between hole dynamics and magnetic 
background in the Cu-0 plane is considered the cen- 
tral issue in high-Tc superconductivity. This interplay is 
quite likely essential not only for the mechanism of hole 
pairing but also for other important features present in 
the cuprates such as the pseudogap phasei and various 
types of phase separated, inhomogeneous, statesi ^i'^i'^i^ 

These problems have been extensively studied by a 
number of analytical and numerical techniques. One 
of the most important approaches adopted is through 
the study of microscopic Hamiltonians, such as the t- 
J model in two dimensions (2D). In particular, numer- 
ical techniques have provided a number of highly reli- 
able results for many of the relevant properties which 
are experimentally measured^ Many of these properties 
are related to the problem of holes and binding of holes 
in an antiferromagnet.^'^'^ These numerical results were 
obtained on relatively small clusters as compared, for 
example, with the antiferromagnetic correlation length 
in the doped region or with the superconducting coher- 
ence length in the superconducting phase. As in these 
numerical studies (see below), many theoretical scenar- 
ios, like the "string"— and "spin bag"— pictures, assume 
the presence of an homogeneous antiferromagnetic back- 
ground. 

However, a number of experimental results indicate the 
presence of various kinds of inhomogeneitiesS*^, in ad- 
dition to the most well-known and controversial ones, 
stripes^, which can not easily be included in numerical 
calculations. Although these inhomogeneities appear in 
some specific compounds, there is another more universal 
type of inhomogeneity, in this case of dynamical nature, 
which are the domain walls (DW) between short range 
AF regions which appear when doping the AF half-filled 
insulator. Inelastic neutron scattering studiesi^ have 
shown that the AF correlation length decreases as dop- 



ing increases. For example, the AF correlation length is 
between 10 and 2QA in the doping range 0.04 < a; < 0.15. 
This finite length of the AF ordered domains is missed in 
most numerical calculations on small clusters but even 
in those cases where the largest distance on the clus- 
ter studied is about or larger than the AF correlation 
length, the numerical technique employed should be able 
to determine these dynamical AF domains and to follow 
the movement of holes within this array of domains. So 
far this kind of analysis has not been performed. The 
closest approach has been done by using density matrix 
renormalization group (DMRG) techniques. ^'^^ However 
in this approach, domain walls appear as static features 
due to the use of open boundary conditions in rectangu- 
lar clusters. 

The effect of domain walls on hole dynamics has been 
discussed specially in the context of stripes.^"' It has been 
realized that stripes are essentially DWs with an almost 
one-dimensional pattern. It is also apparent that holes 
acquire a larger mobility by moving along these DWs 
since there is no cost in magnetic energyiJi It is then 
possible that this same physics of holes moving along 
the boundaries of AF domains be present in all doped 
cuprates, even in those where there is so far no clear 
experimental evidence of stripes just due to the pres- 
ence of finite AF domains. It has been also suggested 
that stripes could be the way in which a strongly corre- 
lated electron system avoids the tendency to macroscopic 
phase separationiM Although the presence of phase sep- 
aration is still controversial in the 2D t-J model from 
numerical studies^ it is likely that the proximity to a 
phase separated state is at the origin of the various kinds 
of inhomogeneities mentioned above. 

The purpose of the present study is then to determine 
the presence of finite AF domains in the 2D t-J model 
and the interrelation between hole movement and do- 
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main walls. A finite temperature conventional quantum 
Monte Carlo technique using the checkerboard decom- 
position is usedtf This technique allows the study of 
larger clusters than the ones accessible by exact diag- 
onalization, and more importantly, although working at 
low hole doping, the number of holes in the system is also 
much larger. Not only most numerical calculations have 
so far involved few holes but also analytical approaches, 
like those mentioned above, are essentially independent 
particle pictures. On the other hand, it has been noticed 
that the system with many holes shows a broad range 
of features, absent in the few holes problem, covered by 
the concept of "topological ordering" In fact, not only 
many holes destroy long range AF order but also the 
experimental result that spin stripes order at lower tem- 
perature than charge stripes is presumably a collective 
effect and not due to single-hole physicsfi^ 

The model studied is the exchange anisotropic t-J 
model: 



^ = - ^ E (4c,. + h.c.) 



(1) 



in standard notation. < = 1, J = 0.35. Square clusters 
with side length 8, 10, 12 and 16 with periodic boundary 
conditions (PBC) and hole doping 0.03 < x < 0.10 are 
considered. Temperature is given in units of t. 

It is possible to argue that an enhanced exchange ZZ 
component appears at an effective level due to inter 
plane exchange coupling. For example, in La2CuO4.11 
(x = 0.14), superconducting below 42 K, there is an 
AF correlation length in the c-direction that couples a 
CU-O2 plane at least with its two (one above, one be- 
low) adjacent planes^ Although smaller in strength than 
other 2D interactions (spin-orbit, Dzialoshinskii-Moriya, 
four-site ring exchange), the inter plane coupling is the 
responsible for a finite temperature transition at half- 
filling, below which the correlation length diverges. In 
any case, we are going to work at very small 7, i.e. close 
to the Ising limit of the exchange interaction, or t-J^ 
model, which is not a realistic case for the cuprates but 
is required to reduce the "minus sign problem" of QMC 
simulations (see below). In addition, at half- filling, the 
isotropic 2D Heisenberg model does not have an AF tran- 
sition at finite temperature (Mermin and Wagner theo- 
rem) but the 2D Ising model has a critical temperature 
of 2.269 J = 0.794 t. In any case, as far as the interaction 
between holes and short-range AF order is concerned, the 
isotropic t-J and the t-J^ models lead to similar resultsn^i 

The "minus sign problem" is a very well-known as- 
pect of QMC simulations and it has been extensively 
analyzedi^^ We give here the outline of the sign calcu- 
lation in order to make our results fully reproducible. 
This calculation consists of two parts. First, a local con- 
tribution which takes into account fermion permutations 
inside each cube. Then a global contribution which is 



due to fermion permutations between different plaque- 
ttes at each time slice. The permutations which appear 
when one hole moves from one side of the cluster to the 
opposite side have not been included. This is expensive 
numerically and it is just a boundary contribution. It has 
been also checked that the influence of this contribution 
is negligible, at least in the range of temperatures con- 
sidered. The results shown below have been obtained by 
averaging at least over six independent runs and at one 
of the lowest temperatures, T = 0.3t, each run consisted 
of 700,000 MC steps after thermalization. Only the re- 
sults corresponding to an average sign larger than 0.01 
were considered in the present study. The relative error 
of the measured quantities is around 0.01, i.e. about the 
size of the symbols used in the plots, except at the lowest 
temperatures reached, where the error bars are two times 
larger. 

In order to identify the AF domains on each Trot- 
ter slice on the square lattice, a variant of the Hoshen- 
Koppelman algorithms^ was implemented. This algo- 
rithm consists in sweeping each slice row by row. At 
each site, its state (z-projection of the spin) has to be 
compared with the one of the previous nearest neighbor 
(NN) site on the same row and with the state of the NN 
site on the previous row. Two NN sites belong to the 
same domain if their spins are antiparallel. Care should 
be taken due to spatial PBC for the last site on each row 
and for the last row on each slice. Once the domains are 
identified the number of sites belonging to it is counted. 
In addition, for all spins belonging to a given domain, 
the number of neighbors belonging to the same domain 
and the number of its neighbors belonging to neighboring 
domains are computed. 

We start with a study of the magnetic ordering in the 
t-Jz limit (7 = 0). In Fig. ^a) it is shown the depen- 
dence with hole doping x of the volume fraction of the 
two largest domains, vi and V2 obtained by dividing the 
number of spins belonging to those domains, Vi and V2, 
by the number of cluster sites. It can be seen that al- 
ready at a high temperature, T « t, there is a large, 
percolating AF domain occupying roughly half the vol- 
ume of the lattice. The second largest domain occupies 
on average roughly one quarter of the lattice, and the rel- 
ative volume of the third largest domain is smaller than 
0.1. One should take into account that for randomly dis- 
tributed ±1 numbers on a lattice (T = 00 limit), undoped 
case, the relative volumes are 0.31, 0.17 and 0.10 respec- 
tively. As the temperature is lowered, at around T 0.7t 
(this value decreases as doping increases) there is a rather 
sharp decrease of the volume of the second largest domain 
with a consequent increase of the largest domain i.e., the 
second largest domain collapses into the first one. This 
change gets smoother as x increases. Below this crossover 
temperature region, the largest ordered domain occupies 
0.95 of the cluster for x — 0.031, decreasing smoothly to 
0.81 for X = 0.12. 

In Fig. nib) the resulting AF order as seen from vi 
and V2 is studied by using a more conventional measure 
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FIG. 1: Results for the 16 x 16 (open symbols), 12 x 12 (full 
symbols), and 10 x 10 (stars) clusters as a function of tem- 
perature, for 7 = and various hole densities as indicated in 
the plot, (a) Volume fraction of the largest AF domain (solid 
lines) and second largest domain (dashed lines), (b) Spin-spin 
correlations at the maximum distance. 



FIG. 2: (a) Volume fraction of the largest AF domain (solid 
lines) and second largest domain (dashed lines) , on the 10 x 10 
cluster and various values of the anisotropy 7 indicated in the 
plot, (b) Same as (a) but for 7 = 0.5 and various cluster 
linear sizes and hole densities, indicated in the plot. 



of magnetic order, i.e. the spin-spin {{SfSj)) correlation 
function at the gain there is 

a steep increase of Smax as T decreases presumably indi- 
cating a finite temperature transition in the bulk limit. 
The temperature of this jump decreases from T « 0.7, 
for X = 0.039 ( 16 X 16 cluster) to T « 0.6, for x = 0.097 
(12 X 12 cluster). 

In Fig. EJa) it is shown the dependence of the vol- 
ume fraction of the two largest domains vi and V2 with 
7. Although up to 7 = 0.25 the results do not depart 
appreciably from those of the Ising limit, at 7 = 0.5 it 
can be observed an important reduction of the volume 
of the largest domain, and hence in the AF ordering of 
the system. This behavior of the 7 dependence is consis- 
tent with the one observed for other quantities as will be 
shown below, and it suggests that the behavior found for 
small 7 at high temperature is going to remain valid at 
lower temperatures (and eventually zero temperature) as 
the isotropic limit is reached. Fig. |2tb) contains a study 
of the dependence of these domain sizes with cluster size 
at 7 = 0.5. The finite size effects are not very impor- 
tant for this anisotropy value and these clusters sizes, 
although it may become more important as the isotropic 
Heisenberg case is reached and for smaller sizesiii No- 
tice however that there is a slightly larger value of the 
doping fraction as the linear size is reduced in the data 
shown in this figure. Hence, if the results were corrected 
at the same doping fraction, finite size effects would be 
somewhat larger. 

Now, the main point of this study is to detect the lo- 
cation of holes in the presence of ordered AF domains 
as long as this doped system is approximately described 
by a t-J model close to the Ising limit. To get some in- 
sight about the microscopic interrelation between holes 



and AF domains it is useful to look at "snapshots" 
of the system generated during the QMC simulations. 
Fig. El^a) shows a picture at a relatively high tempera- 
ture (T = 0.8) where there are still ordered domains with 
the opposite sublattice magnetization (darker shades of 
gray) than the main ordered domain (lightest shade of 
gray). The salient feature in this picture is that holes 
are located preferentially on the boundary of the largest 
AF domain. There are some other situations that can 
be identified in this picture, for example isolated holes 
inside the largest domain with a "string" of overturned 
spins attached to itiAS 

On the other hand. Fig. Otb) corresponds to a typ- 
ical lower temperature behavior. The system presents 
a single ordered domain occupying most of the cluster 
and with just few disordered sites. Some features can be 
identified, such as isolated holes or pairs inside the per- 
colating ordered domain. However, the most important 
feature is that most holes seems to be located in a phase 
separated region. In addition, in this and many other 
snapshots, we observed that most of the spins inside and 
around this region are ferromagnetically aligned, but we 
have not studied this issue systematically. Presumably, 
as AF domains with the opposite magnetization to that 
of the main one shrink with decreasing temperature, they 
would leave behind these highly doped regions forming 
a "mesoscopic" phase separated state or, more properly, 
a charge inhomogeneous state. In this situation holes 
would now gain kinetic energy by moving in a ferromag- 
netic background. 

Of course, it is necessary to determine at a more quan- 
titative level the presence of the features shown in Fig.O 
To this end, we adopted the criterion of looking only at 
the spins and holes belonging to the bulk or to the sur- 
face of the largest ordered AF domain, neglecting the 
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FIG. 3: Snapshot of the 16 x 16 cluster with 16 holes at 
T = 0.8 (top panel) and at T = 0.4 (bottom panel). The 
two lightest shades of gray correspond to spins with the same 
sublattice magnetization, opposite to that of the two darkest 
shades of gray. Holes are indicated by circles. 



contribution from smaller domains. We believe that this 
criterion gives the right behavior of the system. Then, 
for any site belonging to the largest ordered AF domain 
we computed the number of its nearest and next-nearest 
neighbor sites belonging to the same domain. We consid- 
ered a spin in the bulk of the largest AF domain (inner 
spin) if the number of its neighbors belonging to the same 
domain is seven or eight, and we consider it at the bound- 
ary (boundary spin) if that number is between three and 
five. We prefer to leave aside the intermediate case of six 
neighbors and also the case of less than three neighbors, 
in this case to eliminate some sites which are loosely con- 
nected to the main AF domain. In any case, we repeated 
the calculations taking 2-5 neighbors for the boundary 
sites and 6-8 neighbors for the inner sites obtaining es- 
sentially the same results. We are confident then that 
the results shown below are quite robust, independent of 
the details of the classification of spins and holes. 

The average fractions of spins located on the boundary 
(rij) and inside (nf) the largest AF domain, are defined 
as the average number of boundary spins divided by the 
total number of spins in this cluster (Vi) and the aver- 
age number of inner spins divided by Vi respectively. In 
Fig. 2Ja), and nf are shown as a function of tempera- 



FIG. 4: Average fraction of spins located on the boundary 
(dashed lines) and inside (solid lines) the largest AF domain 
as a function of temperature, (a) Same clusters and hole filling 
as in Fig. El 7 = 0. (b) 8 x 8, 6 holes, 7 = (filled triangles), 
0.25 (open triangles) and 0.5 (crosses); 12 x 12, 8 holes, 7 = 
(filled squares), 0.25 (open squares) and 0.5 (pluses). 



ture, at 7 = and other parameters as in Fig. ^ At high 
temperature the fraction of boundary spins is larger than 
the fraction of inner ones but this situation is reversed 
as the temperature is reduced due to the overall order- 
ing of the system. As doping increases, the crossover 
temperature decreases, consistent with the reduction of 
AF order with doping. As it can be seen in Fig. 0fb), 
this crossover temperature is further reduced as the off- 
diagonal exchange term (7) increases. Again, one could 
expect that the region over which there is an important 
fraction of boundary spins extends down to zero temper- 
atures at the isotropic limit. 

To determine the location of holes with respect to these 
AF domains, a similar procedure as above was extended 
to holes. That is, for every hole the number of its nearest 
and next-nearest neighbor spins belonging to the largest 
ordered AF domain was computed. As before, we consid- 
ered a hole located inside the largest AF domain (inner 
hole) if the number of its neighbors belonging to the same 
site is seven or eight, and we consider it at the bound- 
ary (boundary hole) if that number is between three and 
five. The average fraction of boundary holes (nj*) and 
the average fraction of inner holes (n,^) are obtained by 
dividing the number of boundary holes and the number 
of inner holes respectively, by the total number of holes. 

Results of these calculation are depicted in Fig.jS] The 
overall features are similar to those of Fig. 0] i.e. there 
is a dominance of boundary holes at higher temperatures 
changing to inner holes at lower temperatures (Fig.[Sfa)). 
The fact that this crossover temperature at a given dop- 
ing is lower than the corresponding temperature for spins 
suggests that holes have a preference to be located at 
the boundary, a possibility that will be confirmed below. 
Fig. Elb) also shows that this crossover temperature is 
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FIG. 5: Average fraction of holes located on the boundary 
(dashed lines) and inside (solid lines) the largest AF domain 
as a function of temperature, (a) Same clusters and hole 
filling as in Fig. 7 = 0. (b) Same clusters and parameters 
as in Fig.^Ib). 



reduced by increasing 7 again consistent with the results 
of Fig. ^b). Another feature that can be observed is that 
at the lowest temperatures, the reverse behavior occurs, 
i.e., there is an increase (decrease) of the boundary (in- 
ner) hole fractions. It is easy to see that this behavior is 
consistent with the formation of a phase separated hole- 
rich region inside the main AF percolating domain. This 
possibility is also examined below. 

The central problem of the interrelation between holes 
and AF domains is addressed in Fig. El where the ratio 
of the fraction of boundary holes to the fraction of 
boundary spins is shown. This quantity, obtained from 
two quantities computed in QMC simulations, have then 
larger error bars than the ones of those quantities. The 
main result is that the ratio is always noticeably larger 
than one indicating a preference of the holes to move 
along the boundary of AF domains. In Fig. Ela), it is 
shown that increases as T is reduced, and decreases 
as doping increases. In Fig.|SJb), it is also shown that 
decreases as 7 increases. It is out of question to perform 
an extrapolation to 7 = 1 but this behavior suggests that 
one should be cautious about extending the conclusions 
of the present study to the isotropic case. The inset of 
Fig. IHb) shows the hole population of the domain walls, 
Pw , defined as the number of holes on the boundary of the 
largest ordered domain divided by the number of spins 
(i.e., the "length") of this boundary. It turns out that Pw 
is just multiplied by the hole density x and divided 
by vi. To take a reference, on the stripes Pw = 0.5, 
which seems considerably higher than an extrapolation 
to zero temperature of the present data, specially taking 
into account the reduction with increasing 7. 

An elementary argument of minimization of the num- 
ber of broken AF bonds leads to the possibility of pairing 
of holes even when they are constrained to move along the 




0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 



T/t T/t 

FIG. 6: Ratio of the average fraction of holes located on 
the boundary to the average fraction of spins located on the 
boundary as a function of temperature, (a) Same clusters and 
hole filling as in Fig. [2 7 = 0. (b) 8 x 8, 6 holes, 7 = (filled 
triangles), 0.25 (open triangles) and 0.5 (crosses); 12 x 12, 
8 holes, 7 = (filled squares), 0.25 (open squares) and 0.5 
(pluses) . The inset shows hole population of DW for the same 
parameters as in (b), 7 = and 0.5. 



boundary between AF domains. In order to confirm this 
possibility, and to get further indications of the presence 
of a low temperature phase separation, we computed the 
hole-hole correlations C (r) as a function of temperature. 
In Fig.lJfa) these correlations are shown as a function of 
distance and at various temperatures for a fixed value of x 
and 7. As T decreases, the correlations at large distances 
decrease while the correlations at short distances increase 
indicating the onset of hole attraction. Quite interesting, 
at intermediate temperatures, the largest correlation cor- 
responds to the distance of -^2, as was observed in many 
numerical studies for two holes at zero temperature in 
the fully isotropic t-J mode&SiS and this is typical of 
dx^-y^ pairing. At the lowest temperature, the correla- 
tions are enhanced both at near and at far distances and 
suppressed at intermediate distances which is typical of 
a phase separated state. 

For better visualizing the attraction between holes we 
averaged C(r) for r — 1,^/2 and compared it with the 
average over the three largest distances on the cluster. 
Fig. CJb) shows these averages as a function of tempera- 
ture, on the 12 X 12 cluster, 7 = and at various dop- 
ings. It can be seen that the crossover temperature at 
which the near distance correlations start to dominate 
over the long distance ones, i.e. the probability of find- 
ing two holes becomes larger at short distance than at 
long distance, decreases slightly as a function of dop- 
ing (of course, error bars exclude a precise determination 
of this crossover). This crossover temperature could be 
identified with T^^^ as discussed in Ref. IH Fig. [T^c) 
seems also to indicate that the pairing crossover temper- 
ature, for a given cluster size and doping, decreases with 
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FIG. 7: (a) Hole-hole correlations as a function of distance 
on the 12 x 12 cluster, x — 0.069, 7 = 0, at various tempera- 
tures indicated on the plot, (b) Average value of the hole-hole 
correlations at short distance (open symbols) and at large dis- 
tance (filled symbols) on the same cluster, 7 = 0, at various 
hole doping levels indicated in the plot, as a function of tem- 
perature. The arrows indicate approximately the crossover 
temperatures, (c) Same as (b) but for x — 0.069 and various 
values of 7, indicated in the plot. 



increasing 7. One should stress the fact that this onset 
of pairing takes place while there is still a larger frac- 
tion of holes on the DW than inside the largest ordered 
AF domain. For example, from Fig. [7fb), the onset of 
pairing occurs approximately at T = 0.76, 0.70 and 0.66, 
at a; = 0.056, 0.069 and 0.083, respectively. The corre- 
sponding crossover from DW hole fraction to bulk hole 
fraction (Fig.|SJa))takes place at T = 0.69, 0.66 and 0.60 
for the same fillings. 

It is also worth to mention that by looking at hole-hole 
correlations in momentum space we have not found in- 
dications of stripe formation. This hole static structure 
factor is loosely distributed around (0, 0) with not def- 
inite pattern. Presumably, additional mechanisms have 
to be invoked to stabilize a stripe stateii^ The peak of the 
magnetic structure factor is also always at (tt, tt). In addi- 
tion, by examining snapshots like the ones in Fig. |31 and 
by computing a measure of compactness using the data of 
Fig. 0] it turns out that the AF domains are in close con- 
tact between each other and are far from being "round" 
objects as implicitly assumed in Ref. Hence, also in 
this type of inhomogeneity an external agent should be 
present to stabilize the reported static AF islands. 

Let us summarize and discuss the main results of the 
present study. The 2D exchange anisotropic t-J model, 
close to the t-J^ limit, presents at high enough tempera- 
ture finite AF domains which, at least in finite clusters, 
collapse at lower temperatures into a single ordered AF 
state. While there are domains, and hence domain walls, 
holes would tend to preferentially locate along these DW. 
In this case, there are indications of hole pairing starting 



at a relatively high temperature. One could speculate 
that this situation extends down to zero temperature in 
the fully isotropic t-J model as this model is supposed 
to describe, for the appropriate hole doping, the super- 
conducting phase of the cuprates. On the other hand, 
the fully isotropic case may behave as the t-J^ model be- 
haves at zero temperature, at least in finite clusters, and 
in this case it could undergo a crossover to a charge inho- 
mogeneous state. As mentioned in the introductory para- 
graphs, there are many experimental indications of vari- 
ous kinds of inhomogeneities in the cuprates. In any case, 
the nature of the possible charge inhomogeneous state 
detected in the simulations should be further analyzediSi 

A possible consequence of the reduction of the bound- 
ary length as temperature is decreased, is that the holes 
would become increasingly localized. This localization 
could contribute to the increase of in-plane resistivity 
in lightly doped La2-i:Srj:Cu04 and YBa2CuOj, as T 
is reduced reported in recent experimental studies^^*^. 
An indirect support to this scenario comes from the re- 
lation between resistivity and AF correlation length in 
La2_2;Sr2;Cu04. It can be seen (Fig. 3 of Ref. 25) that 
the inverse mobility and the AF correlation length paf 
have an identical behavior as a function of doping. How- 
ever, one should be cautious since the mobility is mea- 
sured at room temperature while paf is the zero tem- 
perature value. Of course, in these compounds, stripes, 
which are a particular kind of DW, could be stabilized. 
In addition, the kinetic energy, computed in our simu- 
lations, is almost constant down to the lowest tempera- 
tures attainable. Perhaps, the system responds to this 
tendency to localization by going into a charge inhomo- 
geneous state. In any case, one could suggest recipes to 
avoid this tendency to localization due to the shrinking 
of AF domains, present in all cuprates. For example, 
DWs could be enforced into the system by introducing 
magnetic impurities outside Cu02 layers in order to pin 
domains with random sublattice magnetizations. The 
resulting enlarged DW length would reduce hole local- 
ization and prevent phase separation. 

Finally, let us compare the present results with the 
ones obtained in some of the most recent and closest 
studies. In the first place, in Ref. 0, DWs are essen- 
tially straight lines, and essentially static, both in their 
analytical treatment and in their DMRG calculations. 
Hence their predictions are of limited applicability in the 
finite temperature real physics of the t-J^ model, as re- 
vealed by the present simulations, where dynamical do- 
mains are considerably non-compact regions. In particu- 
lar, our main result, holes moving essentially along DWs, 
suggests that "transverse" motion of holes is unlikely. 
Notice also that, by looking at the present results on 
a larger "time scale", the AF domains could appear as 
more compact regions separated by broadened disordered 
domain walls. In addition, Ref. is a study at zero 
temperature, where we obtain, with PBC, indications of 
clustering of holes. In this sense, their zero temperature 
DW should be a consequence of additional effects not in- 
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eluded in the t-Jz model. In our case, finite AF domains 
appear as a competition between charge and magnetic 
energy, at finite temperatures. The issue of domain wall 
formation in our model has to be further investigated and 
it would be useful to compare the wall energetics with the 
predictions of Rcf. 15. Another of our most important 
results are the indications of binding of holes. This is a 
delicate issue which has been hotly debated in the con- 
text of stripes li^ii^iiS In Ref. 0it was concluded that the 
presence of a DW is mostly irrelevant to pairing. In stud- 
ies where charge stripes are forced into the system by an 
on-site potential added to an anisotropic t-J model^^, no 
indications of pairing were found: stripes are metallic. In 
these studies, spin and hole correlations along the stripe 
are very much similar to the ones in isolated t-J chains. 



One may suggest then, that the behavior of dynamical 
domain walls studied in the present work is considerable 
different than for stripes where charge is ordered by an 
external mechanism. 
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